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A method is presented to include quantum corrections into the calculation of two-pion corre- 
lations for the case where particles originate from resonance decays. The technique uses classical 
information regarding the space-time points at which resonances are created. By evaluating a sim- 
ple thermal model, the method is compared to semiclassical techniques that assume exponential 
decaying resonances moving along classical trajectories. Significant improvements are noted when 
the resonance widths are broad as compared to the temperature. 
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I. INTRODUCTION 



Analyses of two-boson correlations have provided intangible information regarding the space-time development of 
hadronic reactions [EJ0]- Pions, kaons, and photons have all been exploited for their bosonic nature which results 
in a positive correlation at small relative momentum. Numerous other correlations, involving nucleons |3|-|6| or light 
nuclei Mjq] that are correlated due to the strong or Coulomb interaction as well as identical-particle statistics, have 
also been analyzed and have given further information regarding collision dynamics. Source sizes and time scales have 
been extracted from collisions covering a wide assortment of reactions, from heavy-ion collisions at a few MeV per 

qv . nucleon, where time scales of thousands of fm/c have been determined, to e + e~ — ► jets, where lifetimes of a fraction 

f^*i ' of a fm/c have been observed H . 

The comparison of theoretically predicted correlation functions with experimental results provides an important 
test of the dynamical properties of reaction models. Most models provide semiclassical information about the source 
function S(p,x), the probability of emitting a particle of momentum p from the space-time point x. By convoluting 
the source functions for particles of momenta pi and P2 with the squared relative wave function |^ q (a;i — X2)| 2 , 
one is able to predict the correlation function C(pi,p2). Source functions are usually obtained from semiclassical 
simulations, where the source points are associated with the last point of interaction jlO|]. Particles from resonances 
are usually assumed to be emitted according to an exponential decay law, with the characteristic time usually chosen 
to be independent of the energy of the resonance. Quantum considerations have been explored by Lednicky and 



Progulova 1 1 and by Bertsch, Danielewicz and Herrmann p 

In this study, we pursue two goals. First, we wish to quantify the importance of quantum treatments by comparing 
to semiclassical forms for a simple thermal model. Although the formalism we present is not much different from that 
presented previously in the literature |OJl2 ] , the differences with semiclassical treatments have not been quantitatively 
documented. We find that quantum corrections become important when kinematics constrain the resonances to be 
off shell. Second, we propose an alteration to the methods for extracting correlations from classical simulations to 
better account for quantum effects. We show that a simple modification of the semiclassical treatment can account 
for the quantum corrections by incorporating information regarding the off-shell energy of the decaying resonance. In 
this study we neglect any interaction between the particles aside from the constraints imposed by symmetrization. 

In this paper, we will first briefly review correlation functions for direct sources, which will provide the foundation 



for correlations from resonant sources in the following section. Section IV contains a calculation for a simple Brcit- 
Wigner resonance to demonstrate the importance of quantum treatments. Modified correlation weights are presented 
in Sec. [V] as a means to better calculate correlations from classical simulations when resonant decays are involved. 
Finally, we will compare the various methods mentioned in this paper by considering a simple thermal model in Sec. 

II. REVIEW: CORRELATIONS FROM DIRECT SOURCES 

The two-particle correlation function is usually defined as 



C(pi,P2) 



dN^/(d 3 Pl d 3 P 2 



(dNW/d 3 Pl ) {dNW/d 3 p 2 ) 



(1) 



Before describing the two-particle probability, we first derive an expression for the one-particle probability, which 
also introduces the definition of the source function. One-particle probabilities can be determined by matrix elements 
Tf (x) where / describes the remainder of the system, and x is the point at which the pion has the final interaction 
with the system. Without loss of generality, one can write 



2F dNil) 
The definition of the source function for pions is 



E 



d 4 x T f (x)e ip - X 



S(p, x) = ^ j di5x T f( x + Sx/2)T f (x - 5x/2)e 
which leads to the simple relation 



-ip-Sx 



I 



dNW 
2E p ^^= jd 4 xS(p,x) 



(2) 



(3) 



(4) 



Po=E p 



The source function can be interpreted as the probability per unit space-time for creating a pion of momentum p. 

Since source functions can be extracted from semiclassical simulations or thermal models, it has proven useful to 
also express two-particle probabilities in terms of source functions. The two-particle probability requires a two-particle 
matrix element T\ (x a ,x b )- Assuming independent, or uncorrelated, emission means that the two-particle element 
factorizes, fl2|. 



Tf ] {x a ,X b ) > Tf a (Xa)Tf b {x b ) 

Here, a and b label independent sources. The two-particle probability then becomes 



(5) 



(2E 1 )(2E 2 



dN& 



d 3 p\d 3 p2 



E 

a,b,f a ,fb 



d 4 x a d 4 Xb Tf a (x a )Tf b (x b )—= {exp(ipi • x a + ip2 • x b ) + exp(ip 2 • x a + ip x ■ x b )} 

v2 



(6) 



y^ S d 4 X a d 4 X b Sa{pi,X a )S b (p2,X b ) 



+ d 4 x a d 4 x b S a (p,x a )S b (p,x b )exp[i(p 2 -p{) ■ (x a - x b )] 

^2{Ia(Pl,Pl)h(P2,P2) +I a (Pl,P2)h(P2,Pl)} , 



(7) 
(8) 



a , b 



where p — (jpi +i?2)/2, and the functions I a ,b are defined as 



Ia,b(pi,P2) = I d^x S a , b ( El±Pl 7X ) e ^-P^ 



(9) 



The above allows one to write the correlation function as 



C( P i,p 2 ) = l + 



J2 I a(Pl,P2)h{P2,Pl) 

a , b 

J2 I a(Pl,Pl)h{P2,P2)' 
a,b 



(10) 



This formulation is similar to that originally used by Shuryak |L3] . 

Simulations typically provide a sampling of the on-shell source function. The application of Eq. (0) to simulations 
is made difficult because the source functions are evaluated at po ^ Et pi+P2 \/2 meaning that they require off-shell 
information. The above formalism can be related to simulations through the smoothness approximation [|14|Jl5| , 



S a ( r ' Xa ) ^ b ( O '^J * ^ (P 1 '^)^ 6 iP^i x b) • (11) 

For thermal sources, one can justify a form for the off-shell behavior of S and the smoothness approximation can be 
averted. A truly quantum theory would provide the T matrices that would also allow off-shell evaluation of S, and in 
fact, formalisms have been developed where classical simulations are augmented by converting the point particles into 
wave packets |16| . This also allows one to forego the smoothness approximation, but at the price of inserting an ansatz 
for the quantum behavior that has some peculiarities. We will sidestep the issue of the smoothness approximation in 
this study as we wish to focus on quantum aspects associated with the propagation of off-shell particles. 

When calculating correlation functions from simulations, particles from resonances are usually included in the source 
function by using the space-time points from which the resonances decay. The lifetime of the resonance affects the 
correlation function through the exponential decay that is simulated in the transport model. As will be demonstrated 
later in the paper, the exponential decay law is modified if the dynamics emit resonances with a particular mass or 
range of masses. In this case, the form of the source function becomes nonexponential as will be explained in the next 
two sections. 

III. CORRELATIONS FROM RESONANT SOURCES 

In this section we present a formalism for calculating two-particle correlation functions from resonance decays given 
that one or both of the pions might originate from a resonance. The result will depend on the source function that 
represents the creation of the resonance rather than the source function that represents the points at which the 
final-state pions are created. The evolution and decay of the resonance will be accounted for through the quantum 
propagator of the resonance. The space-time point at which the resonance decayed to produce the final-state pions will 
be treated as an intermediate quantum step in the evolutionary path between the initial creation of the resonance and 
the asymptotic momentum states of the decay products. The expressions derived here thus incorporate an integration 
over all points at which the resonance might have decayed relative to the points at which the resonance is created. 

We will consider three examples, a hypothetical scalar A that decays into two pions, p, which is a vector resonance 
also decaying into two pions, and to, which is also a vector resonance, but decays into three pions. The matrix element 
for creating the pion and a second particle with momentum k through a scalar resonance A is 

T«, fA , k (x)=g [d 4 x A G A {x-x A )T AjA {x A )e lk -\ (12) 



where G A is the Fourier transform of the propagator for the resonance and f A refers to the state of the remainder of 
the source. Here T„ is effectively the T-matrix element for emission of the pion, while T A is the T-matrix element 
that would describe emission of the resonance if the resonance were stable. 

Following the same method as in the previous section, one can use T v to create the pion source function using 
Eq. (|3|) to express T\{x)T A (x') in terms of the source function of the resonance. The resulting expression for the pion 
source function can be used to generate I A , as defined in Eq. (ph, which is all that is needed to calculate correlation 
functions, 



Ia(pi,P2) = g 2 / 7TT d 4 x A exp[i(p 2 -pi) -x A ]S A f 1 2 +k,x A \ G A (pi + k)G* A (p 2 + k) (13) 

G A(P) = , f2 1 TT I 2\ ( 14 ) 

p z - M A + iU A (p z ) 

n A (m 2 ) = M A r A ^^A. (15) 

qrt m 

Here, m 2 = p A , M A and T A are the mass and width of the resonance, respectively. The relative momentum of the 
outgoing pions in the frame of the resonance is q 2 — m 2 /4 — m' 2 r and q^ is the same quantity for an on-shell resonance. 
We emphasize that the interference term can be calculated in Eq. (Oh without reference to the direct source of the 
pions, as the source function of the resonance becomes the required input. 

In order to understand the role of the propagators, we reconsider the case of emitting a pion pair, with momenta p 
and k through a single, scalar resonance with momentum p A = p + k, 



dN 
' d 3 p 



^E p -^ z = I A (p,p) (16) 



dN= f^f^ d i xS A (p + k,x) g - j (17) 

2E p 2E k Ky ! \{ P + kY-M\ + iIi A ({p + kf)\ 2 



A 

d 4 x Sa(pa,x) dm z -lm— 5 — = , 2 . ■ (18) 



— —d 4 x S a (pa,x) dm 2 -Im— —5— 



The spectral function describes the density of states of resonances of invariant masses m, 

11 dn . 

2^ m m 2 _ M 2 _ in A (^) = d^z ■ ( 9 ^ 

Thus, one can see that the source function does not provide any information regarding the mass or width of the 
resonance. However, combined with the spectral function, which derives from the product of propagators, it provides 
the probability of creating the resonance at space time point x with momentum p and with invariant mass m. For the 
direct case considered in the last section, the source function was always evaluated on shell, i.e., the spectral function 
was effectively a delta function. 

One can also perform the same calculation for vector resonances such as the p meson. In that case, where the 
coupling of the pions to the vector meson is 

C bA = iX{nd (l ir)pf l , (20) 

the expression for I p is similar to Eq. ( J13| ) with the source function and propagator accounting for the vector nature 
of the p, 

I p (pi, P 2) = x2 J^r J d4x P ex P[*(P2 ~ Pi) ' *p](Pi - k) a Gf( Pl + k)S p 01 (El±Pl + k} x )j apSfa + k){p2 _ k)s . 

(21) 
The propagator for the vector resonance is 

For the vector case the self-energy scales differently as a function of the resonance mass than in the scalar case, 

n p (m 2 ) = M p r P 4— , (23) 

where the same notation as in Eq. (|l5|) was used. 

As a third example, we consider the propagation of an to meson, which is also a vector resonance, but decays into 
three pions through 

Ant = me^W^TT+dV^TT- . (24) 

In this case the expression for 1^ becomes even more complicated than the p example. 

4i(Pi,Pa) = k 2 / ' 2 ~e~2e / ^ Xu ex P[*(^ 2 ~Pi) ' x »>\ 

xe Q ^c P^k u ^Gf( Pl + k + t)S$, (E±±Pl +k + i }X \ G *7<5 (p2 + k + i) eSflluli , p£Vj«'. (25) 

The expression for the self-energy is also somewhat more complicated. 

/d 3 k d 3 l 
^—— S ((m -E k - E[) 2 - El +l - 2klcos6) m 2 |k x 1| 2 . (26) 

where 9 is the angle between k and 1 and B is an uninteresting constant fixed by the condition that Ii ul {M 2 ) = M^T^, 
Mu and T^ being the mass and the width of u), respectively. After applying the delta function, there remains a fairly 
complicated expression, which was evaluated numerically in our considerations of the lo case later in the paper. 



IV. LIMIT OF A NARROW RESONANCE 



For this discussion we shall, for simplicity, consider the source function for a narrow Breit-Wigner resonance A. 
Then, according to Eqs. (|2|) and (|12|), the probability of emitting a pion pair with momenta k and p is 



(2E k ){2E p )-J^=Y J 9 J<Pxd*XA G A {x - x A )T A j jk (x A )e l ^ + ^ x 



(27) 



= g 2 / d x A d i x S A (p + k,x A )K(p + k,x — x A ). (28) 

Here, K represents the probability of a resonance carrying momentum p + k propagating from x A to x. It may be 
expressed in terms of the propagators, 



K{P ' X)= I$^ e^G* A (p + 5q/2)G A ( P -6 q /2). 



(29) 



In general if is a complicated function. In order to illustrate the quantum nature, we consider the limit of a narrow 
resonance that allows the propagator to be expressed in a simplified form, 

Ga{p) = ^ * . W9 - (30) 

2M A p a - E p + iT/2 

Furthermore, we consider the case where the particles are emitted with equal and opposite momentum. Integrating 
over spatial coordinates gives the probability of the resonance propagating for a time t with off-shellness AE, 

2TM 2 , f o . x r „ . , r , sm(2AEt) . s 

Tl(AE, t) = ± / d 3 x K{p, x) = -e(t)e~ rt \ „ ' , (31) 

7T J IT i\hi 

where AE — E k + E p — M A is the off-shellness. If one integrates over the off-shellness, the expected exponential 
behavior is obtained. 

fd(AE) K(AE, t) = re" rt , (32) 

whereas integrating over t describes the preference for emitting the particle with energy close to on-shell. 

JdtnAE,t)= l - AE2 r f {T/2)2 . (33) 

The oscillating term sm(2AEt), which is responsible for preferentially emitting resonances with small AE, also 
governs the distribution of emission times. Classical simulations, which are typically based on Monte Carlo algorithms, 
cannot easily incorporate regions with negative probabilities as suggested by this form. The mean propagation time 
in a transport simulation could be altered to match the mean time of the quantum propagator, 

jdttn(AE,t) r/2 

W ' JdtTZ(AE,t) AE 2 + (T/2) 2 ' { ' 

However, the second moment for the time would not correspond to that expected for an exponential decay with the 
same average time. 

1 3r*-4A^ 

8 (AE 2 + T 2 /4) 2 W V ' 

Not only is this form inconsistent with exponential decay, (t 2 ) can even become negative for resonances far off 
shell. This illustrates, on a formal level, the need for performing the quantum corrections described in this paper. 
Quantitative comparisons follow in the next section. 



V. CORRELATION WEIGHTS 

Simulations of heavy-ion collisions usually provide the creation points of pions along with their outgoing momenta. 
Neglecting other interactions besides symmetrization, correlation weights for pions of momenta p\ and p 2 originating 
from space-time points x a and Xb are usually determined by calculating the average symmetrization weight of all pairs 
satisfying the imposed binning or acceptance pC| . 

Y J Jd i X a d i X b S a (pi 1 X a )S b (p2 1 X b )w a (pi,p 2 )w b {p2,Pl) 

C(pi,p 2 ) = l + - ^ r , 4 -T A 7T7 T7T7 s (36) 

2_, J d^Xad^Xb b a (pi,X a )b b (p2,Xb) 
a, b 

= 1 + (w a (p 1 ,p 2 )w b (p 2 ,pi)). (37) 

Inspection of Eqs. (0), (flQ), and (11) reveals the weight for direct sources, 



gl Pl+P2 x \ 

Wd{pi,P2) = exp[z(pi -p 2 ) -x d ] 2 — '-— (38) 

Wd (Pi>P2) = exp[i(pi -p 2 ) ■ id], (39) 

where the semiclassical form assumes the smoothness approximation |Ij,n5[ . 

If the particle originates from the decay of a scalar resonance A, the weight takes on a different form as can be seen 
from inspecting Eq. (113), 

, x r 7 X 1 S A (S1±E2- + k X A ) ( Pl + fc) 2 M\ - t n A [( Pl + fc) 2 ] 

w A (pi,P2) =exp»(pi ^2-ia g 7T ,. w — — TT^ jTTf TT77 — Ira' ( 40 ) 

sU(pi + fc, au) (pa + fc) - ^| - «n j4 [(P2 + fc) 2 ] 

Here the resonance was created at x A and decayed into pions of momenta p\ and fc. The space-time coordinate of the 
decay does not enter the weight as all decay points have been considered. 

If the decay is of a vector resonance such as a p meson, the weights are slightly different, 

, x ,., s 1 s p { Pl , P 2, fc) (pi + fc) 2 - m 2 - »n p [( Pl + fc) 2 ] 

«,,(*,*) = «*[•(* -p 2 ) • ^^—-—-^-^j^^^-—, (41) 



S p (pi,P2,k,x p ) = (pi -k) a S p ap I 1 - +fc,x p j (p 2 -fc)' 5 . 



(42) 



In the derivation of iSp we have assumed the two pions involved in the resonance decay to have equal mass. 

Any sort of resonance can be included in this manner, including resonances that decay into three or more bodies. 
One such example is lo, which decays into three pions, one of each species. Labeling the momenta of the two pions, 
whose symmetrization we ignore, as fc and I, the weights turn out to be 

f v r. f v , S u (p UP 2,k,l) (Pl + fc + If - Ml - iUvKp! + fc + Q 2 ] 

Wu,{pi,P2) = expi pi -p 2 ) -Xu, \-=- — t- 2 —2 rpj-fT n21 , (43) 

Su(Pi,Pi,k,l) (p 2 + k + iy - M£ - i\l u \{p2 + k + I) 2 \ 

S u (pi,P2,k,l,x») = e% £ Pi*"^ {^ El + k + l 1 x}j e\, v>v p%'k»'t '. (44) 

Thus, weights can be used to calculate correlation functions for the decay of any resonance in a rather straightforward 
manner. The formalism coherently accounts for all points at which the resonance may have decayed, but requires 
information regarding the points at which the resonances were created as well as information about the accompanying 
particles in the decay. The only difficulty comes in assigning the ratios of the source functions, i.e., the smoothness 
problem. 

For a thermal source, S A (p,x) ~ poe~ po < T in the scalar case and S^Jp,x) ~ poe~ po ' T (—g a p + PaPp/p 2 ) in the 
vector case. In the thermal cases the Boltzmann factor cancels out of the ratios used to calculate weights. Thus 
for the thermal example, the weights become a product of four factors, a phase arising from the points at which the 
resonance R (= A, p,u>) is created, a ratio of energies, a spin factor, and a ratio of propagator denominators. 

(t/i), x rf , i (^ + k ') ■ n Xfl(pi, P2 ) (Pi + k') 2 - M 2 R - l n R [( Pl + fc') 2 ] 

w R Pi,P2 = exp«(pi -p 2 ) • Xr] , , ,,s t w -ijto slf TT77 TTT^T' ( 45 ) 

(Pi + fc')'« XMPi,Pi) (P2 + fc') 2 - M j?. - ?r T?i(P2 + fc') 2 ] 



where n refers to the frame of the thermal source and k' is equal to either k for a two-body decay, or to k + I for a 
three-body decay. 

The spin factors xr depend on the sort of resonance being considered. 

Xa(pi,P2) = 1 (46) 

Xp(pi,Pa) = (Pi - k) a (~g af 3+p a pp/p 2 ) (pa - k) p , (47) 

X*(pi,J*) = ^ riWl* {-g a p+p a p /p 2 ) e^' p^W, (48) 

where p = (pi +p 2 )/2 + k' . 

However, one cannot easily apply this technique to nonthermal sources because the off-shell behavior of the source 
function is not always known. This problem also confronts calculations with direct sources, and forces one to either 
invoke the smoothness approximation, or assume some form for the off-shell behavior of the source function. 

In analogy to the smoothness approximation one might use the thermal weight w R assuming a reference frame n^ 
or simply neglect the ratio of energies in Eq. (Eq). In fact, the ratio X.r(pi,P2)/x,r(pi,Pi) can also be neglected to a 
reasonable approximation as the structure of the correlation function derives largely from the last factor in Eq. J4q), 
the ratio of propagator denominators. 

VI. COMPARISON WITH SEMICLASSICAL MODELS 



In this section we compare calculations of the quantum type described in Sec. Ill with semiclassical calculations 
where the resonance is assumed to propagate classically and decay according to an exponential form exp(— t/r). For 
the purposes of this comparison we choose to model two simplified systems, one of decaying p resonances and a second 
of decaying w's. For each case resonances are produced and decayed with a Monte Carlo procedure according to a 
thermal distribution characterized by a temperature of 150 MeV. The mass of p is 770 MeV and the width is chosen 
to be 150 MeV, while the mass of w is 783 MeV and the width is 8.4 MeV. 

To model the uncorrelated emission of pion pairs, we thermally create particles of momentum k (and I for the w) 
by Monte Carlo, then add to the particles the weight 

w . = E k{+ E l} + E i XRJPtiPi) , 4g ., 

E k {Ei}Ek \{ pi + k 'f-M 2 R + in[{p i + k') 2 ]f 

where i — 1,2 and the braces indicate terms that appear for u) only. This weight accounts for the spectral function 
of the resonances as described in Eq. (|19|) and for the spin factors \R as m Eqs. (fig)- (fig). We note that the weights 
W1.2 are merely used to generate resonances and their products, and are not related to the correlation weights. If 
these weights were included through a keep/reject prescription, they would not need to appear in any of the following 
expressions. 

Once the pions are statistically generated, one simply calculates the average weights described previously to generate 
the correlation functions. 



C(p llP2 ) = l + K 



(Ew 1 a wf^(p 1 , P2 )Ew^w^ b (p 2 , Pl y 

a b 



V 



EWT EW 3 6 



(50) 



For our comparison wr is either the weight for p, w pi or the one for u>, w u . Note that a and b label individual 
resonances. 

In the semiclassical descriptions, the weights are determined by calculating the expectation value 

w^\ Pl , P2 ,x R ) = (e^-^ x ) (51) 

= exp[i(pi - p 2 ) ■ xr] / (^(x - x R ) S 3 [x - x_r - v R (t - t R )] exp[z(pi - p 2 ) ■ (x - x R )] 

x exp[-(* - t R )/(j R T)]Q(t - t R ) (52) 

IRJ 

rv ^ i m R /TR , 

= exp[i(pi - pa) • x R ] — r, (53) 

mR/TR + ip R ■ (pi - pa) 



where 7^ is the Lorentz factor due to the motion of the resonance. Here, -u/ sc ' assumes an exponential form for the 
pion emission, which is characterized by a lifetime r. The same form for exponential decays was developed by Padula 
and Gyulassy fll7| ]. 

Based on one's perspective, one might choose any of several prescriptions for the energy dependence of the lifetime 
T(m). We investigate three possibilities: (1) The lifetime is chosen such that m/r = II(r7i 2 ). This choice would be 
motivated by the form of the propagator. (2) The lifetime is chosen to correspond to the average emission time as 



described in Sec. IV, except that the relativistic generalization of Eq. (|34|) is used r = 2m 5(m 2 — M| — ill) -1 . (3) 
A fixed lifetime 1/T is used. 

If the resonance distribution function in coordinate space is independent from that in momentum space, the inter- 
ference term in the correlation function factorizes into a term stemming from the space-time extent of the resonance 
source itself and one that arises from the decay process. 

C( Pl , P2 ) - 1 = »{(exp[t(pi -pa) • (x A - x B )]) (C'(p 1 ,p 2 ) - 1)} (54) 

Here, xa and xb refer to points at which the resonances are created, and (exp[i(pi — p 2 ) ■ (xa — Xb)]) represents the 
weighted average using the product of the source functions as the weight. 

(exp[i(pi - p 2 ) ■ (xa - x B )}) = \Jr(pi - pz)\ (55) 

Y, f d 4 x A S A [(pi + Pi) /2 + k A , x A ] cxp[i(pi -p 2 ) ■ x A ] 

Jr(Pi -Vl) = — ^ f , 4 FT-J7 i wo 1 ; i ( 56 ) 

2^ J d 4 x A SA[{Pi +P2)/2 + kA,x A \ 

A 

The reduced correlation function, C'(pi,p 2 ) — 1, is similar to the average of the weights (w a Wb) from Sec. [VJ, only with 
the factors of exp[i(pi — p 2 ) ■ xr] removed from the weights in Eq. (Ea). Since C" contains all the relevant information 
about the decay, we will focus on the reduced correlation function for our comparisons. 

Figures |l| and g display the reduced correlation function for the cases where (pi + p 2 )/2 = 200 MeV/c and 800 
MeV/c, respectively. The upper panel of each figure displays results for the case where pi is parallel to p 2 while the 
lower panel displays the results for the case where the two momenta are perpendicular. All three semiclassical results 
exhibit significant deviations from the quantum calculations. 

Figures || and display the same information but for a thermal source of u> mesons. In this case, the semiclassical 
treatments are significantly more accurate. This was expected as the width of the u) is much less than the temperature, 



which allows the spectral function of the u) to be sampled evenly. As described in Sec. IV, the distribution of decay 
times become exponential when evenly averaged over all masses. 

The overall width of the correlation functions in the upper panels of Figs. |3| and ^ can be understood by noting 
that the correlation's width should be determined by the condition AEt = 1, where AE s=s VRq. The width of the 
correlation function for the lower panels, where pi _L p 2 , is more complicated since it is more sensitive to the spatial 
movement of the resonance while it decays. The correlation functions for the p in Figs, w and ^ are even more 
complicated since they extend to large relative momenta where AE approaches q. Given the complicated kinematics 
involved, it is not surprising that the result is sensitive to the exact form for the semiclassical treatment. 

In the case of the p, none of the three semiclassical prescriptions for r(m) provides a consistently good approximation 
to the quantum result for all scenarios shown in Figs, [fl andu. They all exhibit significant deviations from the quantum 
result for high relative momentum of the pion pair, above 200 MeV/c. In the case of the to, where only small relative 
momenta, q < 200 MeV/c, are relevant, prescription (1) seems to best reproduce the quantum result, as can be seen 
in Figs. I and § 

It should be emphasized that the correlation from the non-zero extent of the resonance source function has been 
factored out in this calculation. The deviations of the semiclassical results occur for relative momenta of a hun- 
dred MeV/c or more. In a heavy-ion collision the factor (exp[i(pi — p 2 ) ■ (xa ~ Xb)]) hi Eq. ( |55| ) would tend to zero 
for relative momenta much greater than 50 MeV/c due to the large spatial sizes of the emitting regions. Thus, the 
form of C becomes irrelevant for higher relative momenta unless the source sizes are small, e.g., pp collisions. 

VII. CONCLUSION 

Our findings imply that the semiclassical treatments work quite well for the larger sources considered with heavy-ion 
collisions. However, for smaller sources, especially when the resonance widths are comparable to the temperature and 
resonances are produced far off-shell, the semiclassical treatments significantly deviate from the quantum result. This 
can be linked to the failure of the usual exponential decay law when off-shell resonances are involved as was shown in 



Sec. IV 



In Sec. [V| we proposed a method to correctly calculate correlation functions from semiclassical models, exploiting 
the source function of the resonance and its creation point in space-time as opposed to the creation points of the 
final-state pions. This method could be easily applied to generate correlation functions from the event histories of 
simulations. When using direct pions, the creation points of the final-state pions provide all the necessary information 
for creating correlation functions. By considering the creation points of the resonances that decayed into the final- 
state pions, one is able to coherently account for all space-time points at which the resonance might have decayed by 
modifying the prescription for generating correlation weights as described in Sec. M. 

It should be emphasized that the quantum formalism discussed here becomes important only for sources that 
themselves are quantum-mechanical in nature. If the source is large and the product of the momentum and spatial 
uncertainties are large, Ap Ax ^> h, the behavior of the correlation function is dominated by the term e tq '( XA ~ XB ' , 
which is determined by the points at which the resonance is created. Hence C plays little role in the correlation 
function of the p resonance described in Fig. YU when the overall source size is many Fermi as is the case for heavy- 
ion collisions. Quantum considerations in resonant decays could play an important role when considering the decay 
of small sources that push the limits of the uncertainty principle. However, such sources are also accompanied by 
questions regarding the quantum nature of the source functions responsible for initial creation of the resonances, i.e., 
the smoothness approximation might not be justified. For such problems, unless one knows the off-shell behavior of 
the source functions as is the case for a thermal model, the quantum treatments presented here address only half the 
problem. 
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FIG. 1. Reduced correlation function for a source of p mesons at a temperature of 150 MeV. The average momenta of the 
two pions is fixed at 200 MeV/c. By factoring out the space-time dependence of the p, the manifestation of the p lifetime is 
singled out. The exact quantum treatment is shown to differ from the three semiclassical treatments that are described in the 
text. The failure of the semiclassical descriptions owes itself to the fact that the p width is sufficiently large for the thermal 
source to effectively emphasize a specific region of off-shellness. 
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FIG. 3. Same as Fig. |l[ but using a thermal source of lu's at a temperature of 150 MeV. In this case, since the width of the 



u> is much less than the temperature, semiclassical treatments work remarkably well. 
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FIG. 4. Same as Fig. H, except that w's have average momentum P — 800 MeV/c. 
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